
function [wage] = w_mat(det_earn, i,j,k, dt )
wage=exp(det_earn(1) + det_earn(2).*i + det_earn(3).*i.*i + det_earn(4).*i.*i.*i + det_earn(5).*j + det_earn(6).*j.*j + k).*dt; 
end
